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Abstract 

LUNASKA (Lunar UHE Neutrino Astrophysics with the Square Kilometre Array) is an ongoing project conducting 
lunar Cherenkov observations in order to develop techniques for detecting neutrinos with the next generation of radio 
telescopes. Our current observing campaign is with the 64-metre Parkes radio telescope, using a multibeam receiver with 
300 MHz of bandwidth from 1.2-1.5 GHz. Here we provide an overview of the various factors that must be considered 
in the signal processing for such an experiment. We also briefly describe the flux limits which we expect to set with our 
current observations, including a directional limit for Centaurus A. 

Keywords: UHE neutrino detection, coherent radio emission, lunar Cherenkov technique, UHE neutrino flux limits, 
detectors-telescopes 
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1. Introduction 

The lunar Cherenkov technique is a method for detect- 
ing ultra-high energy (UHE) particles which makes use of 
the regolith of the Moon as a target volume, giving it a 
larger potential aperture than any other current method. 
A UHE particle shower in this material will produce a co- 
herent radio pulse through the Askaryan effect PQ, which 
may be observed by terrestrial radio telescopes, as pro- 
posed by Dagkesamanskii and Zheleznykh [2]. Active ex- 
periments of this type include NuMoon g], RESUN g], 
and our own LUNASKA (Lunar UHE Neutrino Astro- 
physics using the Square Kilometre Array) [5]. 

The long-term goal of the LUNASKA project, as the 
name suggests, is to develop strategies for and eventually 
utilise the Square Kilometre Array: a large radio telescope 
currently being planned for future construction in either 
Australia or South Africa Our observations thus far 
have made use of the Australia Telescope Compact Array 
(ATCA) _S| - and, now, the Parkes radio telescope. Here 
we focus on the signal processing requirements of our cur- 
rent observations, the understanding of which is essential 
for lowering the energy threshold at which UHE particles 
may be detected. 
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2. Observations with the Parkes Radio Telescope 

The Parkes radio telescope is a single-dish telescope 
with a diameter of 64 metres, located in New South Wales, 
Australia. We are conducting 200 hours of lunar Cherenkov 
observations with it during 2010, using its 20cm multibeam 
receiver. This receiver has 300 MHz of bandwidth from 
1.2-1.5 GHz, and allows us to place multiple pointings or 
'beams' on the sky simultaneously, with dual linear polar- 
isations. Each beam has a full width at half maximum 
(FWHM) size of around 14 arcminutcs. 

A radio pulse from a lunar Cherenkov event is expected 
to come from the limb of the Moon, with radial polari- 
sation, which guides our decision on where to point the 
beams. In addition, the system noise in our observations 
is dominated by the thermal noise of the Moon, so we must 
compromise between pointing more on the Moon (less sen- 
sitive), or slightly off the Moon (more sensitive; but ob- 
serving a smaller area). We have explored several pointing 
strategies, but have spent most observing time in the con- 
figuration shown in figure [I] This configuration allows us 
to keep two beams slightly off the Moon, with their linear 
polarisations aligned such that the power of a single event 
may be detected mostly on a single polarisation channel. 

The sensitivity of the lunar Cherenkov technique to 
UHE particles is highly directionally dependent [7], which 
makes it possible to target observations at a particular 
potential source, to maximise sensitivity to it. In our ob- 
servations, we have targeted the radio galaxy Centaurus A, 
which is associated with an excess of UHE cosmic ray ar- 
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Figure 1: Pointing strategy with Parkes radio telescope. Circles rep- 
resent the positions of beams with respect to the Moon, and crosses 
within them indicate the alignment of the linear polarisations. The 
off-moon beam is for anticoincidence filtering - see section [3. 4| 



rival directions [5]. We note that the total power of the 
neutrino flux from Centaurus A is limited by gamma-ray 
observations [5]; however, the contribution of the flux in 
our energy range to the total power is small, so it is less 
strongly constrained. 

3. Signal processing 

We use the standard front-end signal path for the 20 cm 
multibeam receiver at Parkes, with minor modifications 



(see section 3.1|. The receiver amplifies the signal, and 
mixes it with a sine-wave local oscillator at 1555 MHz to 
downconvert it from the 1.2-1.5 GHz radio band to 55-355 
MHz. 

At this point, the signal enters our purpose-built digital 
backend, the 'Bedlam' Neutrino Transient Processor. It is 
converted from analog to digital with 8 bits of precision, 
at a rate of 1.024 Gsamples/s. Then, it is passed through 
a finite impulse response (FIR) filter, with 64 coefficients 
or 'taps' (see section 3.2). The samples, both before and 
after filtration, are stored in a running buffer with a length 
of up to 8k samples. A threshold test is applied both to 



these and to interpolated values (see section 3.3) and, if 
the appropriate anticoincidence criteria arc met (see see- 
the contents of all buffers are written to disk. 



3.4) 



tion 

While the buffers are being stored, the system is not pro- 
cessing incoming data: the fraction of this 'dead time' was 
generally held to around 5%. 



Further processing (see sections 3.5 and 3.6) is applied 



to the stored buffers in post-analysis. The real-time pro- 
cessing is not required to be exhaustive: it only needs 
to ensure that any potentially significant lunar Cherenkov 
event is stored, so that it can be identified later. 

3.1. Attenuation 

We tested the response of the telescope by transmitting 
a sine- wave from a small antenna on the telescope surface. 
The received sine-wave was observed to be clipped, with 



its peaks truncated. This is similar to the receiver sat- 
uration problem encountered in a similar experiment us- 
ing another telescope, the Very Large Array, by Jaeger et 
al [2]. In our case, the problem was traced to a particu- 
lar amplifier which was being driven outside of its design 
voltage range, and failing to maintain linear amplification 
of the signal. This was solved by attenuating the signal 
earlier in the signal chain, bringing the amplifier back into 
its design range, and using our own amplifiers to restore 
the signal afterwards. Subsequent repetition of the earlier 
test confirmed that this had the desired effect. 

3.2. Dedispersion 

It has been recognised since the first lunar Cherenkov 
observations by Hankins et al [10] that a broadband co- 
herent pulse from the Moon will be dispersed as it passes 
through the ionosphere, causing the low-frequency com- 
ponents to be delayed with respect to the high-frequency 
components, and thus reducing the peak amplitude of the 
pulse. In our previous experiments with the ATCA, we 
used an analog de-dispersion filter to compensate for this 
effect. In our current observations, we use a digital finite 
impulse response (FIR) filter implemented in our digital 
backend. This has the advantage that the filter can be 
easily reconfigured to apply a different amount of dedis- 
persion, so we can adjust it during each observing session 
to match changes in ionospheric conditions and the slant 
depth of our sight-line to the Moon. A typical filter con- 
figuration is shown in figure [2j and its performance char- 
acterised in figure [3] 

The parameter of the ionosphere that determines the 
magnitude of the dispersion is its total electron content 
(TEC), for which we obtain hourly figures from the Aus- 
tralian Bureau of Meteorology [TT], who produce iono- 
spheric models based on ionosonde radar measurements. 

3.3. Interpolation 

If the sampling rate of the digital signal were infinite, 
then any pulse would always be sampled at its peak, allow- 
ing a threshold test to detect its full amplitude. However, 
the finite sampling rate of any real system means that 
there will not generally be a sample taken exactly at the 
peak of the pulse, and so the full amplitude of the pulse 
will not be detected. The loss of sensitivity due to this 
effect is described in [5]. 

Since the input signal is band-limited, and sampled 
at over the Nyquist frequency, the Shannon-Nyquist sam- 
pling theorem tells us that our samples contain full in- 
formation about the original analog signal, including any 
missed peaks, and the Whittaker-Shannon interpolation 
formula allows us to obtain missing intermediate values 
by convolving the signal with a sine function. This al- 
lows the full amplitude of a pulse to be recovered, and is 
computationally cheaper than directly increasing the sam- 
pling rate. However, it is still somewhat computationally 
intensive, and becomes more so as more intermediate val- 
ues are required, so our real-time hardware reconstructs 
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Figure 2: Impulse response of de-dispersion filter compensating for 
7 ns of dispersion, which is a typical maximum value during one of 
our observing sessions. This filter is equivalent to correlating with an 
expected pulse shape the same as this impulse response. Note the 
separation of high-frequency components (left) and low-frequency 
components (right). 

only one intermediate value between each pair of real sam- 
ples (thus effectively doubling the sampling rate) - which 
is sufficient to ensure that any significant real pulse will be 
stored - and more exhaustive interpolation is performed in 
post-processing. 

3.4- Anticoincidence 

A major concern for observations of this type is that 
there exists radio-frequency interference (RFI) that resem- 
bles nanosecond-scale pulses of the sort that would be ex- 
pected from a lunar Cherenkov event. Such pulses can be 
detected from the surroundings of the telescope via the 
far sidelobes of its beam. However, since we are using a 
multibeam receiver, these pulses are distinguished by ap- 
pearing on all beams simultaneously. By contrast, a real 
lunar Cherenkov event should appear only on the beam 
pointing at the part of the Moon where it occurred. 

We therefore apply an anticoincidence filter, excluding 
potential events which appear in more than one beam. A 
simple such filter is applied in real time, excluding the ma- 
jority of pulse-like RFI, and a more sophisticated version 
is applied in post-analysis. The off-moon beam, having a 
lower system temperature, is more sensitive, and particu- 
larly effective at identifying RFI pulses. 

3.5. Phase search 

Recent simulations of Cherenkov pulses in the time- 
domain show that their shape is bipolar, with separate 
positive and negative peaks [IS]. In the frequency do- 
main, this corresponds to the pulse being entirely imag- 
inary, with a phase of ±7r/2. However, the process by 
which a radio telescope converts a signal from radio fre- 
quency (RF) to a lower intermediate frequency (IF) will 
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Figure 3: Performance of filter shown in figure [2] measured as devi- 
ation from the ideal delay at each radio frequency. The band con- 
taining the signal is 1.2-1.5 GHz. 

change this phase randomly, depending on the phase off- 
set between the arriving pulse and the local oscillator (LO) 
signal. We implement a search through the possible shapes 
of the original pulse, for different possible LO phases. 

3.6. Band shaping 

Optimum detection of a known pulse against a back- 
ground of noise requires 'pre- whitening' the noise spec- 
trum, then correlating with the expected pulse shape |13j . 
In the frequency domain, this is equivalent to weighting 
different frequencies according to the expected signal-to- 
noise ratio at each of them. 

The thermal noise from the Moon follows a Rayleigh- 
Jeans spectrum, with spectral power proportional to v 2 , 
which is the same as the spectrum of a fully-coherent 
Cherenkov pulse. Since this source dominates the total 
noise present in the signal, the signal-to-noise ratio is ap- 
proximately constant across all frequencies, and the entire 
band should be weighted equally. However, it may be de- 
sirable to weight lower frequencies more heavily to opti- 
mise for the detection of partially-incoherent Cherenkov 
pulses, for which the spectrum cuts off at higher frequen- 
cies. 

4. Sensitivity 

Using the Monte Carlo simulation described in [14] , we 
calculate the sensitivity of our experiment, and the lim- 
its we expect to set to a flux of UHE neutrinos (assuming 
no detections). For an isotropic flux, our projected limit 
is shown in figure [4j along with limits from other exper- 
iments in the same energy range. For a directional flux 
from Centaurus A, our projected limit is shown in figure 
\E\ Note that Centaurus A has extended structure on a 
scale comparable to the size of the region on the sky from 
which we are sensitive to UHE particles, so the limit to a 
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Figure 4: Differential limits (90% confidence) on the isotropic all- 
flavour neutrino flux set by lunar Cherenkov observations (thick 
lines) and other experiments (thin lines). The curve for Parkes repre- 
sents the limit we expect to set with our current observations, assum- 
ing no detections in 200 hours. Other lunar Cherenkov observations 
include RESUN A [3], NuMoon [3], and our own observations with 
the ATCA [5]. Other experiments include RICE [15], ANITA (T6| . 
and the Pierre Auger Observatory (PAO) | 17 | (for upgoing neutrinos 
only). A 1:1:1 flavour ratio is assumed. 

flux from its outer structure (i.e., its giant lobes) will be 
slightly different. 

The version of the simulation software used here docs 
not incorporate any variation of the gradient of the lunar 
surface over the scale of a particle shower. The effect of 
surface roughness on this scale will generally be to reduce 
the aperture of the experiment at low energies, and to 
increase it at high energies; see appendix B of [5] for more 
discussion of this effect. 

5. Conclusion 

We are currently conducting our most sensitive lunar 
Cherenkov experiment to date. We are deliberately tar- 
geting Centaurus A, and expect that these observations 
will set the lowest limit on its neutrino flux in the energy 
range to which we are sensitive. We are refining the signal- 
processing techniques that will be used in future experi- 
ments with more sensitive telescopes, eventually including 
the SKA. 
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Figure 5: Differential limits (90% confidence) on an all-flavour neu- 
trino flux from Centaurus A. As in figure [4] the curve for Parkes is 
for the projected limit from our current observations. Also shown 
are the limit from our previous observations with the ATCA |18j . 
a limit derived 1181 from the exposure of RICE 1151 . and a prelimi- 
nary integrated (not differential) limit set by the PAO |17| . Other 
experiments in figure [4] have minimal exposure to Centaurus A. 
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